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Abstract — Cognitive  radio  networks  (CRNs)  have  the  potential  to  utilize  spectrum  efficiently  and  are  positioned  to  be  the  core 
technology  for  the  next-generation  multi-hop  wireless  networks.  An  important  problem  for  such  networks  is  its  capacity.  We  study 
this  problem  for  CRNs  in  the  SINR  (signal-to-interference-and-noise-ratio)  model,  which  is  considered  to  be  a  better  characterization 
of  interference  (but  also  more  difficult  to  analyze)  than  disk  graph  model.  The  main  difficulties  of  this  problem  are  two-fold.  First,  SINR 
is  a  non-convex  function  of  transmission  powers;  an  optimization  problem  in  the  SINR  model  is  usually  a  non-convex  program  and 
NP-hard  in  general.  Second,  in  the  SINR  model,  scheduling  feasibility  and  the  maximum  allowed  flow  rate  on  each  link  are  determined 
by  SINR  at  the  physical  layer.  To  maximize  capacity,  it  is  essential  to  follow  a  cross-layer  approach;  but  joint  optimization  at  physical 
(power  control),  link  (scheduling),  and  network  (flow  routing)  layers  with  the  SINR  function  is  inherently  difficult.  In  this  paper,  we  give 
a  mathematical  characterization  of  the  joint  relationship  among  these  layers.  We  devise  a  solution  procedure  that  provides  a  (1  -  e) 
optimal  solution  to  this  complex  problem,  where  s  is  the  required  accuracy.  Our  theoretical  result  offers  a  performance  benchmark 
for  any  other  algorithms  developed  for  practical  implementation.  Using  numerical  results,  we  demonstrate  the  efficacy  of  the  solution 
procedure  and  offer  quantitative  understanding  on  the  interaction  of  power  control,  scheduling,  and  flow  routing  in  a  CRN. 
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1  Introduction 

OGNITIVE  radio  networks  (CRNs)  have  great  po¬ 
tential  to  improve  spectrum  efficiency  and  will  be 
fhe  core  fechnology  for  fhe  nexf-generafion  mulfi-hop 
wireless  nefworks.  Wifhin  such  a  nefwork,  each  node  is 
equipped  wifh  a  cognitive  radio  (CR)  for  wireless  com¬ 
munications,  which  is  capable  of  sensing  fhe  available 
frequency  bands  (i.e.,  fhose  bands  fhaf  are  currenfly  nof 
used  by  primary  users)  and  reconfiguring  RF  fo  swifch 
fo  newly-selecfed  frequency  bands. 

From  wireless  nefworking  perspective,  CR  offers  a 
whole  new  sef  of  research  problems  in  algorifhm  design 
and  protocol  implemenfafion.  To  appreciafe  such  oppor- 
funify  we  compare  CR  wifh  a  closely  related  wireless 
fechnology  called  multi-channel  multi-radio  (MC-MR)  [2], 
[20],  [21].  Firsf,  MC-MR  employs  fradifional  hardware- 
based  radio  fechnology  and  fhus  each  radio  can  only 
operafe  on  a  single  channel  af  a  time.  In  confrasf,  fhe 
radio  technology  in  CR  is  software-based;  a  soff  radio  is 
capable  of  swifching  frequency  bands  on  a  per-packef 
basis.  As  a  resulf,  fhe  number  of  concurrenf  frequency 
bands  ufilized  by  a  CR  is  f5rpically  much  larger  fhan  by 
MC-MR.  Second,  a  common  assumption  for  MC-MR  is 


•  Y.  Shi  and  Y.T.  Hou  are  with  the  Bradley  Department  of  Electrical 
and  Computer  Engineering,  Virginia  Polytechnic  Institute  and  State 
University,  Blacksburg,  VA,  24061.  E-mail:  {yshi,thou}@vt.edu. 

•  S.  Kompella  is  with  Information  Technology  Division,  U.S. 
Naval  Research  Laboratory,  Washington,  DC,  20375.  E-mail: 
sastry.kompella@nrl.navy.mil. 

•  H.D.  Sherali  is  with  the  Grado  Department  of  Industrial  and  Systems  En¬ 
gineering,  Virginia  Polytechnic  Institute  and  State  University,  Blacksburg, 
VA,  24061.  E-mail:  hanifs@vt.edu. 

Manuscript  received  June  10,  2009;  revised  Jan.  3,  2010  and  March  27 ,  2010; 

accepted  July  1,  2010. 


fhaf  fhere  is  a  sef  of  "common  channels"  available  af  ev¬ 
ery  node  in  fhe  nefwork.  Such  assumpfion  is  hardly  frue 
for  CRNs,  where  each  node  may  have  a  differenf  sef  of 
available  frequency  bands.  These  imporfanf  differences 
befween  MC-MR  and  CR  suggesf  fhaf  algorifhm  design 
for  CRNs  is  subsfanfially  more  complex  fhan  fhaf  for 
MC-MR  nefworks.  In  some  sense,  an  MC-MR  nefwork 
can  be  considered  as  a  special  case  of  a  CRN. 

The  capacify  of  a  CRN  is  usually  considered  a  key 
problem  in  fundamenfal  undersfanding.  In  fhis  paper, 
we  aim  fo  sfudy  fhis  problem  in  fhe  SINR  (signal-fo- 
inferference-and-noise-rafio)  model.  In  fhis  model,  con¬ 
currenf  fransmissions  are  allowed  and  interference  (due 
fo  fransmissions  by  non-infended  fransmiffer)  is  freafed 
as  noise.  A  fransmission  is  successful  if  and  only  if  SINR 
af  fhe  receiver  is  greater  fhan  or  equal  fo  a  fhreshold. 
Moreover,  fhe  achieved  fransmission  capacify  is  also  a 
function  of  SINR  (via  Shannon's  formula).  The  currenf 
undersfanding  is  fhaf  fhe  SINR  model  is  better  fhan  fhe 
so-called  "disk  graph  model"  (or  "protocol  model"  [15]) 
for  inferference  characferizafion. 

Alfhough  fhe  SINR  model  is  fhoughf  fo  be  more 
realisfic  fhan  fhe  profocol  model,  a  number  of  problems 
arise  when  frying  fo  carry  ouf  mafhemafical  analysis 
in  if.  Firsf,  SINR  af  a  receiver  nof  only  depends  on 
fhe  fransmission  power  af  fhe  corresponding  fransmiffer, 
buf  also  depends  on  fhe  fransmission  powers  af  ofher 
fransmiffers.  Mafhemafically,  SINR  is  a  non-convex  func- 
fion  of  multiple  variables.  Many  optimization  problems 
in  fhe  SINR  model  are  non-convex  problems  and  NP- 
hard.  Second,  since  bofh  scheduling  feasibilify  and  fhe 
maximum  allowed  flow  rate  on  each  link  are  determined 
by  SINR,  an  optimal  solution  fo  maximize  capacify  musf 
be  developed  wifh  joinf  consideration  of  network,  link. 
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and  physical  layers.  Due  to  these  difficulties,  theoretical 
results  on  CRNs  in  the  SINR  model  remain  limited. 

In  this  paper,  we  investigate  the  network  capacity 
problem  for  mulfi-hop  CRNs  in  fhe  SINR  model.  For 
a  given  insfance  where  each  node  has  access  fo  a  sef 
of  available  bands  (likely  heferogeneous),  we  sfudy  fhe 
network  capacify  problem  via  a  cross-layer  opfimizafion 
approach.  In  parficular,  we  consider  how  fo  maximize 
fhe  rafes  of  a  sef  of  user  communication  sessions,  wifh 
joinf  considerafion  af  physical  layer  (via  power  confrol), 
link  layer  (via  frequency  band  scheduling),  and  nef- 
work  layer  (via  flow  routing).  We  give  a  mafhemafical 
characferizafion  of  fhese  layers  and  formulafe  a  mixed 
infeger  nonlinear  program  (MINLP)  problem.  For  fhis 
opfimizafion  problem,  we  firsf  identify  core  opfimizafion 
variables  and  core  opfimizafion  space  based  on  fhe  phys¬ 
ical  significance  of  fhe  variables.  We  devise  an  algorifhm 
based  on  fhe  branch-and-bound  framework  fo  obfain  a 
(1  —  e)  optimal  solution. 

Alfhough  branch-and-bound  framework  is  sfandard, 
many  componenfs  wifhin  fhis  framework  are  nof  spec¬ 
ified.  We  design  fhe  following  problem-specific  com¬ 
ponenfs.  (1)  We  propose  a  reformulafion-linearizafion 
fechnique  (RLT)  fo  develop  a  fighf  linear  relaxafion  so  as 
fo  obfain  a  fighf  upper  bound.  (2)  For  fhe  lower  bound, 
we  design  a  local  search  algorifhm  by  analyzing  and 
removing  infeasibilify  in  fhe  linear  relaxafion  solution. 
(3)  For  problem  partitioning,  we  propose  fo  choose  a 
parfifion  variable  based  on  ifs  physical  significance.  Wifh 
these  problem-specific  designs,  the  overall  solution  can 
find  a  (1  —  e)  optimal  solution  much  faster  than  brute- 
force  exhaustive  search. 

Our  theoretical  result  offers  a  solution  for  a  given 
network  instance  and  can  serve  as  a  performance  bench¬ 
mark  for  a  multi-hop  CRN.  Such  a  result  is  not  yet 
available  in  the  literature.  If  the  available  bands  at 
each  node  in  the  network  varies  on  a  relatively  larger 
time  scale  than  the  execution  of  our  solution,  then  our 
solution  can  be  used  to  compute  a  (1  —  e)  optimal  result 
for  each  time  period  during  which  the  available  bands 
at  each  node  is  static.  If  the  available  bands  at  each  node 
in  the  network  varies  on  a  relatively  smaller  time  scale, 
then  the  performance  benchmark  can  be  obtained  offline 
by  our  solution.  One  can  use  our  result  as  a  reference 
performance  benchmark  to  measure  the  quality  of  some 
other  proposed  algorithms  (possibly  heuristic,  despite 
distributed)  that  are  developed  for  actual  deployment. 

The  remainder  of  this  paper  is  organized  as  follows. 
In  Section  2,  we  review  related  work  on  cross-layer 
optimization.  Section  3  gives  a  mathematical  character¬ 
ization  of  power  control,  scheduling,  and  routing  in 
the  SINR  model  for  multi-hop  CRNs.  In  Section  4,  we 
reformulate  the  optimization  problem  and  obtain  a  clean 
and  compact  formulation.  Section  5  analyzes  the  core 
optimization  space  and  describes  the  main  algorithm  to 
obtain  (1  —  e)  optimal  solution.  In  Section  6,  we  develop 
tight  upper  and  lower  bounds  inside  the  main  algorithm. 
We  discuss  how  to  interpret  and  apply  our  theoretical 
result  in  Section  7.  Section  8  presents  numerical  results 
and  Section  9  concludes  this  paper. 


2  Related  Work 

There  have  been  active  research  efforts  on  cross-layer 
optimization  for  wireless  networks.  Many  of  these  efforts 
are  based  on  the  disk  graph  model  (see,  e.g.,  [2],  [18], 
[20],  [24]  for  wireless  networks  with  traditional  radios 
and  [26],  [28],  [33],  [34]  for  CRNs  in  particular).  Under 
this  model,  the  impact  of  interference  from  neighboring 
nodes  is  solely  determined  by  whether  or  not  a  node 
falls  within  the  interference  range  of  other  transmit¬ 
ting  nodes.  The  controversy  surrounding  (or  arguments 
against)  the  protocol  model  is  that  a  binary  decision  of 
whether  or  not  interference  exists  (based  on  interference 
range)  does  not  accurately  capture  physical  layer  charac¬ 
teristics.  As  a  result,  the  accuracy  (and  validity)  of  results 
developed  under  protocol  model  remains  unclear. 

In  contrast,  the  SINR  model  is  widely  regarded  as  a 
better  model  for  interference  characterization.  Although 
such  model  is  preferred,  there  are  many  difficulties  in 
carrying  out  analysis  with  this  model  due  to  the  com¬ 
putational  complexity  SINR  involves,  particularly  when 
we  study  cross-layer  optimization  in  a  multi-hop  envi¬ 
ronment.  As  a  result,  many  previous  efforts  were  done 
on  single-hop  networks,  e.g.,  [3],  [10],  [13],  [14],  [17].  For 
multi-hop  networks,  various  simplifications  have  been 
employed.  For  example,  in  [4],  the  authors  assumed 
S5mchronized  power  control,  where  transmission  power 
at  each  node  in  the  network  is  adjustable  but  is  "s5mchro- 
nized"  (identical).  Needless  to  say,  such  s5mchronization 
in  power  control  cannot  offer  optimal  network  perfor¬ 
mance.  There  are  also  some  efforts  studying  cross-layer 
problems  involving  two  layers  instead  of  three  layers 
(physical,  link,  and  network,  as  we  do  in  this  paper). 
For  example,  in  [6],  Bhatia  and  Kodialam  optimized 
power  control  and  routing,  but  assumed  some  frequency 
hopping  mechanism  is  in  place  for  scheduling,  which 
helps  simplify  joint  consideration  of  scheduling.  In  [9], 
Elbatt  and  Ephremides  optimized  joint  power  control 
and  scheduling,  but  assumed  routing  was  given  a  priori. 
In  [27],  Shu  and  Krunz  studied  how  to  maximize  the 
total  rate  on  all  links  in  a  CRN,  with  the  consideration 
of  power  control  and  channel  assignment.  Eor  cross-layer 
optimization  in  the  SINR  model  involving  three  layers 
(physical,  link,  and  network),  nearly  all  existing  efforts 
(e.g.,  [7],  [8])  followed  a  "layer-decoupled"  approach  to 
simplify  analysis.  Under  such  an  approach,  the  final  so¬ 
lution  is  obtained  by  determining  algorithm  /  mechanism 
for  one  layer  at  a  time  and  then  piecing  up  them  together 
instead  of  solving  a  joint  optimization  problem.  Due  to 
de-coupling  in  the  solution  procedure,  these  approaches 
are  heuristic  at  best  and  cannot  offer  any  performance 
guarantee. 

On  another  line  of  research,  various  efforts  have  been 
made  to  study  asymptotic  behavior  (or  scaling  laws)  of 
wireless  networks  (see,  e.g.,  [1],  [5],  [15],  [16],  [19],  [21], 
[22],  [29],  [30],  [31],  [32]).  These  efforts  differ  from  ours 
in  this  paper,  which  focuses  on  designing  optimal  cross¬ 
layer  algorithms  for  a  finite-sized  network. 

3  Mathematical  Models 

We  consider  a  multi-hop  network  with  a  set  of  CR  nodes 
Af.  Each  node  i  G  Af  senses  its  environment  and  finds  a 
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set  of  available  frequency  bands  A4i  for  fhe  given  fime 
insfance  (i.e.,  fhose  bands  fhaf  are  currenfly  nof  used 
by  primary  users),  which  may  nof  be  fhe  same  as  fhe 
available  frequency  bands  af  ofher  nodes.  We  assume 
fhaf  fhe  bandwidfh  of  each  frequency  band  (channel)  is 
W.  Denofe  A4  fhe  union  of  all  frequency  bands  among 
all  fhe  nodes  in  fhe  network,  i.e.,  A4  =  IJieA^  Denofe 
A4ij  =  which  is  fhe  sef  of  frequency  bands 

fhaf  is  common  befween  nodes  i  and  j  and  fhus  can 
be  used  for  communicafion  befween  fhese  two  nodes. 
In  the  rest  of  this  section,  we  present  mathematical 
characterization  of  each  layer  and  obtain  a  formulation 
for  the  capacity  problem  in  this  study. 


3.1  Power  Control,  Scheduling,  and  Their  Relation¬ 
ship  in  the  SINR  Model 


Power  control  on  each  transmitting  node  at  the  physical 
layer  affects  SINR  at  a  receiving  node.  These  SINR  values 
in  turn  will  affect  scheduling  decision  at  link  layer.  That 
is,  if  a  node  is  scheduled  to  receive,  then  its  SINR  must 
be  at  least  a  (minimum  requirement).  Therefore,  power 
control  and  scheduling  are  tightly  coupled  via  SINR  and 
cannot  be  modeled  separately.  In  this  section,  we  formu¬ 
late  power  control,  scheduling,  and  their  relationship  in 
the  SINR  model. 

Scheduling  for  transmission  at  each  node  in  the  net¬ 
work  can  be  done  either  in  frequency  domain  or  time 
domain.  In  this  paper,  we  consider  scheduling  in  fre¬ 
quency  domain  in  tbe  form  of  assigning  frequency  bands 
(charmels).  Since  time  domain  based  formulation  is  sim¬ 
ilar  to  that  for  frequency  domain,  our  approach  can  be 
extended  to  time  domain  based  formulation  and  will 
yield  similar  results. 

To  maximize  the  capacity,  there  may  still  be  concurrent 
transmissions  within  the  same  channel  (and  thus  inter¬ 
ference).  Denote  scheduling  variables  a;™  as  follows. 


1  If  node  i  transmits  data  to  node  j  on  band  m, 
0  otherwise. 


Due  to  interference,  a  node  i  can  use  a  band  m  for 
transmitting  to  only  a  single  node  j  or  receiving  from 
only  a  single  node  k.  That  is, 

m&Mk  meMj 

^  xTi+  1  {i£N',me  Mi)  .  (1) 


For  power  control,  we  assume  that  the  transmission 
power  at  a  node  can  be  tuned  to  a  finite  number  of 
levels  between  0  and  Pmax-  To  model  this  discrete  power 
control,  we  introduce  an  integer  parameter  Q  that  repre¬ 
sents  the  total  number  of  power  levels  to  which  a  trans¬ 
mitter  can  be  adjusted,  i.e.,  transmission  power  can  be 
0,  ^Pmax,  ^Tlnax,  •  •  •  jT’max-  Denote  gj”  G  {0,  1, 2, •  •  •  , Q} 
the  integer  levels  for  transmission  power.  Clearly,  when 
node  i  does  not  transmit  data  to  node  j  on  band  m, 
is  0.  Thus,  power  control  and  scheduling  is  coupled  with 
each  other  via  the  following  relationship. 


e[l,g]  Ifa:™  =  l, 

=  0  otherwise. 


{i,j£Af,i^j,meMij)  (2) 


Consider  a  transmission  from  node  i  to  node  j  on  band 
m.  When  there  is  interference  from  concurrent  transmis¬ 
sions  on  the  same  band,  SINR  at  node  j,  denoted  as  s™, 
is 


{i,j  j,m£  M,j)  ,  (3) 

where  rj  is  the  ambient  Gaussian  noise  density  and  gij 
is  the  propagation  gain  from  node  i  to  node  j. 

Note  that  in  theory,  for  any  small  SINR,  the  corre¬ 
sponding  capacity  is  still  positive  (by  Shannon's  capacity 
formula).  But  in  practice,  if  SINR  is  too  small,  then  the 
achieved  capacity  will  also  be  very  small.  In  this  case, 
such  a  weak  link  will  not  be  very  useful  to  carry  traffic 
flow.  Thus,  we  may  use  a  threshold  to  remove  such  weak 
links  from  considerations.  In  this  regard,  we  introduce 
a  threshold  for  SINR,  i.e.,  a  transmission  from  node  i 
to  node  j  on  band  m  is  considered  successful  if  and 
only  if  s™  >  a.  Then  we  have  the  following  coupling 
relationship  for  scheduling  (xf))  and  SINR  (s™). 

x™  =  1  s“  >  a  (i, .  (4) 

3.2  Routing  and  Link  Capacity 

We  assume  there  is  a  set  of  £  active  user  communication 
(unicast)  sessions  in  the  network.  Denote  s{l)  and  d{l)  the 
source  and  destination  nodes  of  session  I  G  £  and  r{l) 
the  minimum  rate  requirement  (in  b/s)  of  session  L  In 
our  study,  we  aim  to  maximize  a  common  scaling  factor 
K  for  all  session  rates.  That  is,  we  aim  to  determine  the 
maximum  K  such  that  a  rate  of  K-r{l)  can  be  transmitted 
from  s{l)  to  d{l)  for  each  session  /  G  £  in  the  network. 

To  route  each  data  flow  from  its  source  node  to  its 
corresponding  destination  node,  multi-hop  relaying  may 
be  necessary,  due  to  limited  transmission  power  at  each 
node.  Further,  for  optimality  and  flexibility,  it  is  desirable 
to  allow  flow  splitting  and  multi-path  routing.  This  is 
because  a  single-path  flow  routing  tor  a  session  is  overly 
restrictive  and  is  unlikely  to  offer  optimal  solution. 

Mathematically,  this  can  be  modeled  as  follows.  De¬ 
note  fij{l)  the  data  rate  on  link  z  — >  j  that  is  attributed 
to  session  1.  If  node  i  is  the  source  node  of  session  /,  i.e., 
i  =  s{l),  then 

£  /,,(0  =  r(0N  {l€£,t  =  sil))  .  (5) 

If  node  z  is  an  intermediate  relay  node  for  flow  attributed 
to  session  I,  i.e.,  z  f  s{l)  and  z  f  d{l),  then 

E  /h(o=  E 

j^J\f  ,j^i,s{l)  kGAf,k^i,d{l) 

{le£,ieJ\f,ifs{l),d{l)).  (6) 

If  node  z  is  the  destination  node  of  session  /,  i.e.,  z  =  d{l), 
then 

fUl)=r{l)K  (lGAz  =  ci(0)  .  (7) 
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In  addition  to  the  above  flow  balance  equations  at  each 
node  i  G  J\f  for  session  I  e  C,  the  aggregated  flow  rates 
on  each  radio  link  cannot  exceed  this  link's  capacity.  For 
a  link  f  — >  j,  we  have 

E  /h(o<  E  w\og^{i  +  s^) 

lec  meMij 

(i,  j  j,  Mij  ^  0)  .  (8) 

This  constraint  shows  the  coupling  relationship  between 
flow  routing  and  SINK. 

3.3  The  Capacity  Problem 

In  this  paper,  we  study  a  capacity  problem  for  mulfi-hop 
CRNs.  For  capacify  problem,  fhe  simplesf  objecfive  is  fhe 
sum  of  fhroughpuf  achieved  by  all  sessions.  However, 
such  an  objecfive  may  lead  fo  biased  rafe  allocafion 
among  sessions  [23].  Anofher  objecfive  is  maxmin,  i.e., 
maximizing  fhe  minimum  fhroughpuf  among  all  ses¬ 
sions.  Maxmin  has  been  used  in  a  number  of  works  (e.g., 
[15],  [21])  fo  ensure  fairness.  In  fhis  work,  we  consider 
how  fo  maximize  a  common  scaling  factor  K  for  all 
sessions  under  some  given  minimum  rafe  requiremenfs 
r{l),  i.e.,  whaf  is  fhe  maximum  facfor  K  such  fhaf  a  rafe 
of  K  ■  r{l)  can  be  fransmiffed  from  s{l)  fo  d{l)  for  each 
session  I  G  C  in  fhe  network.  Nofe  fhaf  fhis  objecfive  is 
a  more  general  form  of  maxmin  in  fhe  sense  fhaf  when 
r{l)  =  1  for  each  session  I  G  C,  fhis  objecfive  becomes 
fhe  maxmin  fhroughpuf  objecfive.  On  fhe  ofher  hand, 
when  r{l)  ^  1,  i.e.,  fhe  minimum  required  rafe  for  each 
session  I  may  be  differenf  from  session  fo  session,  fhe 
objecfive  funcfion  becomes  fo  maximize  each  session's 
rafe  proportional  fo  ifs  minimum  rafe  requiremenf. 

Puffing  fogefher  fhe  objecfive  and  all  the  constraints 
for  power  control,  scheduling,  and  flow  routing,  we  have 
the  following  formulation. 

Max  K 

s.t.  constraints  (l)-(8) 

x'S  G  {0, 1},  G  {0, 1,  2,  •  •  •  ,  Q},  tT,s7]>0 

{i,j  GjV,iy^j,mGMij) 

K,fijil)>0  {lG£,i,jGJ^,i^j,i^d{l),j^s{l),Mij^9)  . 

Note  that  we  assume  Mi,  the  set  of  current  available 
bands  at  node  i,  is  given  for  a  particular  time  instance. 
The  solution  to  the  above  optimization  problem  will 
offer  the  best  performance  for  this  given  instance.  Since 
Mi  may  change  over  time,  the  optimal  solution  may  also 
change  over  time. 

4  Reformulation 

The  formulation  in  Section  3.3  is  still  in  a  "raw"  form  and 
is  the  first  step  in  setting  up  our  cross-layer  optimization 
problem.  Much  work  needs  to  be  done  to  put  it  into 
a  more  "clean"  and  compact  form  that  is  amenable 
to  mathematical  operation.  In  this  section,  we  analyze 
each  constraint  carefully  and  perform  this  necessary  and 
important  reformulation  step. 


•  The  constraint  described  in  (2)  is  not  suitable  for 
mathematical  programming.  We  reformulate  it  with 
the  following  linear  constraint. 

X™  <  g]"  <  Qx™  {i,3GM,i^],mGM,j)  .  (9) 

It  is  easy  to  verify  that  this  constraint  is  equivalent 
to  (2). 

A  closer  look  at  (9),  in  conjunction  with  (3)  and  (4), 
suggests  that  (9)  can  be  further  simplified  to  remove 
redundancy.  That  is,  we  find  that  (3)  and  (4)  yield 
a;™  <  9™  and  thus  the  left-hand-side  (LHS)  of  (9) 
can  be  removed.  To  show  this,  we  consider  the  two 
cases  for  x^.  For  the  case  when  =  0,  x™  <  g™ 
holds  by  the  definition  of  g™.  For  the  case  when 
x™  =  1,  by  (4),  s™  must  be  positive.  Then  by  (3),  gj” 
must  also  be  positive.  By  the  definition  of  g™,  g™  >  1 
when  it  is  positive.  Therefore,  g™  >  x™.  So  (9)  (or 
(2))  can  be  replaced  by  the  following  constraint. 

q7^<QxTj  {i,jGAf,i^j,mGM^j).  (10) 

•  Constraint  (3)  is  in  the  form  of  a  fraction.  In  a  math¬ 
ematical  programming,  product  is  a  better  form.  We 
can  re-write  it  as 


a  —P 

yij  Q  max 


rjW^ 


l^kGAr,kiCi  l^hGAf,hiCk  ykj  Q  a^max 

_ 9^jq^j _ 

riWQ  I 

Pmax  l^kejG,kiCi  l^hejG,h^k  ykj^kh 

{i,j  G  Af,i  j,m  G  M,j)  . 


This  is  equivalent  to 

U/D  meMk  meMh 

^  m  I  \  '  \  '  m  m  _  m  _  n 

p  .  9kjqkh^ij  9ijqij  —  0 

keAf,kiCi  h£M,h=jCk 

{i,j  G  Af,i  ^  j,m  G  Mij)  .  (11) 

Note  that  in  (11),  g[!j,  and  are  variables  while  all 
other  symbols  are  constants.  Thus,  we  have  a  double 
sum  of  nonlinear  terms  q'^h^Yj  in  (11).  To  reduce  the 
number  of  nonlinear  terms,  denote 

mGMh 

=  E  {k  G  Af ,m  G  Mk)  ■  (12) 

Then  (11)  can  be  re-written  as 


yWQ 

-frnay 


E  9kytTsTy-9^JqY]  =  0 


kGAf  ,k^i 


{i,jGJ\f,i^j,mGMtj),  (13) 


which  now  only  involves  a  single  sum  of  nonlinear 
terms  t'^sYj-  Note  that  by  introducing  new  variables 
t™,  we  decrease  the  number  of  nonlinear  terms  from 
0(|M|4  •  |M|)  in  (11)  to  0(|At|3  •  |M|)  in  (13). 

•  Similar  to  (2),  the  constraint  described  in  (4)  is  not 
suitable  for  mathematical  programming.  We  now 
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Max  K 

.r’"  <  1 

l^keAf,k^i  -^ki  l^j(=J\r,j7^i  ij  - 

(i  e  W,  m  e  Aii) 

qTj  -Q^7j<0 

{i,j  e  Af,i  ^  j,m  &  Mij) 

™  _  n 

(i  e  W,  m  e  Mi) 

Q  1  ^  4-m  ^  _  n 

Pmax  ii  l^k£jq',k^iy>=Pk  ^ij  ^ 

{i,j  e  W,  i  ^  j,m  &  Mij) 

ax™  -  s™  <  0 
— 

{i,j  e  Af,i  ^  j,m  &  Mij) 

{lec,i  =  s{l)) 

-  J2k€'Mtk^i,d(l)  /fedO  =  0 

(le  s{i),d{i)) 

/o(0  -  +  ^[^)  <  0 

(i,j  3,  Mij  7^  0) 

xY]  e  {0,1}, e  {0,l,2,---,Q},t™,s”}  >  0 

(i,j  G  M,i  ^  j,m  G  Mij) 

KJijil)  >  0 

(1  G  C,i,j  G  J\f,i  j,i  d{l),j  i=-  s(l),Mij  i=-  0) 

Fig.  1.  Problem  formulation. 


show  that  (4)  can  be  eliminated  once  we  have  (10), 
(13)  and  the  following  new  constraint. 

s™  >  aa;™  (14) 

We  now  verify  fhaf  (10),  (13),  and  (14)  can  lead  fo 
fhe  following  two  relationships  described  in  (4),  i.e., 

=  1  s™  >  a  (15) 

s™  >  a  x™  =  1  {i,jGAf,i^j,m€Mij)  (16) 

Firsf  we  consider  fhe  case  when  =  1.  By  (14),  we 
have  s™  >  a,  i.e.,  (15)  holds.  Now  consider  fhe  case 
when  s™  >  a.  By  (13),  since  f™'s  are  non-negafive 
variables  and  ri,W,Q,Pmax,gkj,gij  are  all  positive 
consfanfs,  g™  musf  be  posifive.  Then  by  (10), 
musf  be  posifive  when  is  posifive.  Since  x^  can 
only  fake  0  or  1  by  ifs  definifion,  x™  musf  be  1.  Thaf 
is,  (16)  holds. 

•  Finally,  we  can  prove  fhaf  (7)  is  redundanf  once  we 
have  (5)  and  (6).  Thus,  we  will  leave  (5)  and  (6)  in 
fhe  formulafion  and  remove  consfrainf  (7). 

Wifh  fhese  careful  reformulations,  we  now  have  a 
cleaner  and  more  compacf  problem  formulafion,  which 
is  shown  in  Fig.  1. 

5  Solution  Overview 

For  fhe  optimization  problem  in  Fig.  1,  Q,rj,W,a,  Pmax, 
gij,  and  r{l)  are  consfanfs  and  and 

fij{l)  are  opfimizafion  variables.  This  formulafion  is  a 
mixed  infeger  non-linear  program  (MINLP),  which  is 
NP-hard  in  general  [11]  and  cannof  be  solved  by  CPLEX. 
In  Section  5.1,  we  firsf  analyze  fhe  infricafe  relafionship 
among  fhe  variables  and  identify  fhe  core  variables 
among  all  fhe  variables.  We  show  fhaf  fhe  dependenf 
variables  can  be  derived  once  fhese  core  variables  are 
fixed.  We  call  fhe  opfimizafion  space  for  fhe  core  vari¬ 
ables  as  fhe  core  opfimizafion  space.  In  Section  5.2,  we 
presenf  fhe  main  algorifhm  on  how  fo  defermine  an 
opfimal  solufion  in  fhe  core  opfimizafion  space. 


5.1  Core  Variables  and  Their  Optimization  Space 

For  the  complex  MINLP  problem,  its  variables  include 
and  K,  which  collectively  con¬ 
tribute  a  seemingly  huge  optimization  space.  However, 
a  closer  investigation  of  these  variables  show  that  they 
are  inter-dependent.  In  particular,  we  find  that  a;™  and 
variables  are  "core"  variables  and  other  variables  can 
all  be  derived  by  core  variables.  As  a  result,  we  can  focus 
our  study  on  an  optimization  space  by  the  core  variables, 
which  is  a  much  smaller  space. 

We  now  show  how  to  derive  dependent  variables 
from  core  variables.  For  f]"  and  s™  variables,  they  can 
be  derived  by  (12)  and  (13),  respectively,  once  x^  and 
are  given.  For  fij{l)  variables  and  K  variable,  their 
optimal  values  can  be  determined  by  a  linear  program 
(LP).  That  is,  once  we  have  the  values  for  all  a;™  and 
g’”  variables  and  we  have  computed  the  values  for 
all  and  s™  variables  by  (12)  and  (13),  respectively, 
the  optimization  problem  reduces  to  a  network  flow 
problem,  which  is  an  LP. 

With  this  new  understanding  of  optimization  vari¬ 
ables,  we  can  now  focus  our  efforts  on  the  core  variables 
x^  and  q^. 

5.2  Main  Algorithm 

In  this  section,  we  describe  our  main  algorithm,  which  is 
based  on  the  branch-and -bound  framework  [25].  Those 
readers  who  are  familiar  with  branch-and-bound  can 
skip  this  section  and  go  to  Section  6,  where  we  de¬ 
sign  problem  specific  algorithms  for  each  component  of 
branch-and-bound. 

Under  branch-and-bound,  we  aim  to  provide  a  (1  —  e) 
optimal  solution,  where  e  is  a  small  positive  constant 
reflecting  our  desired  accuracy  in  the  final  solution.  In 
case  we  set  e  =  0,  an  optimal  solution  can  be  obtained. 

We  can  start  by  developing  upper  and  lower  bounds 
{UB  and  LB)  for  the  objective  function  (see  Fig.  2(a) 
for  an  example).  Branch-and-bound  requires  to  develop 
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Original  Problem  1 


Problem  2  Problem  3 


Problem  2  Problem  4  Problem  5 


(b)  Iteration  2. 


(c)  Iteration  3. 


(a)  Iteration  1. 

Fig.  2.  Illustration  of  branch-and-bound. 

Initialization: 

1.  Let  the  initial  best  solution  =  0  and  the  initial  lower 
bound  LB  =  —oo. 

2.  Determine  initial  value  set  for  each  core  variable. 

3.  Initialize  the  problem  list  to  include  the  original  problem, 
denoted  this  problem  as  problem  1. 

4.  Obtain  an  upper  bound  UBi  for  problem  1. 

Iteration: 

5.  Select  problem  z  that  has  the  maximum  UBz  among  all 
problems  in  the  problem  list. 

6.  Update  upper  bound  UB  =  UBz. 

7.  Find  a  feasible  solution  t/jz  with  a  lower  bound  LBz. 

8.  If  {LBz  >  LB)  { 

9.  Update  tpe  =  i’z  and  LB  =  LBz- 

10.  If  LB  >  (1  —  e)UB,  we  stop  with  a  (1  —  e)  optimal 
solution  i/jj. 

11.  Otherwise,  remove  all  problems  z'  with  LB  > 

(1  —  e)UB^i  from  the  problem  list.  } 

12.  Build  two  new  problems  zi  and  22  from  problem  2. 

13.  Remove  problem  2  from  the  problem  list. 

14.  Obtain  UBzi  and  UBz2  for  problems  21  and  22. 

15.  If  LB  <  {l  —  e)UBzi,  add  problem  21  into  the  problem  list. 
If  LB  <{l  —  e)U  Bz2f  add  problem  22  into  the  problem  list. 

16.  Go  to  the  next  iteration. 


Fig.  3.  Main  algorithm. 

problem  specific  algorithms  for  finding  upper  and  lower 
bounds,  which  will  be  described  in  Section  6.  After  we 
obtain  the  two  bounds,  we  compare  the  gap  between 
them.  If  the  upper  and  lower  bounds  are  close  to  each 
other,  i.e.,  LB  >  {1  —  e)UB,  then  the  feasible  solution 
corresponding  to  the  current  lower  bound  LB  is  (1  —  e) 
optimal  and  we  are  done. 

Otherwise,  we  need  to  further  narrow  the  gap  between 
UB  and  LB.  To  do  this,  branch-and-bound  partitions 
the  original  problem  1  into  two  new  problems  2  and  3 
(see  Fig.  2(b)).  This  is  accomplished  by  choosing  an 
appropriate  core  variable  x™  or  <7™  and  dividing  its  value 
set  into  two  smaller  sets.  The  choice  of  specific  core 
variables  is  important  as  its  affects  complexity.  We  will 
show  how  to  do  this  in  Section  6. 

After  partitioning,  the  core  optimization  space  is  di¬ 
vided  into  two  sub-spaces  for  problems  2  and  3,  re¬ 
spectively.  We  again  obtain  upper  bounds  UB2  and  UB3 
and  lower  bounds  LB2  and  LB^  for  problems  2  and  3, 
respectively.  Since  the  optimization  space  of  problems  2 
and  3  are  both  smaller  than  that  of  problem  1,  we 
can  have  tighter  upper  bounds,  i.e.,  max{{7i?2,  C/Sa}  < 
UBi,  which  in  turn  yield  better  lower  bounds  with 
max{Li?2,  >  LBi  for  our  maximization  problem 

(see  Fig.  2(b)).  Then  the  upper  bound  of  the  original 


problem  is  updated  as  UB  =  max{[/i32,  CtSa}  and  the 
lower  bound  is  updated  as  LB  =  m.ax{LB2,  LB^}.  As 
a  result,  we  now  have  a  smaller  gap  between  UB  and 
LB.  Then  we  either  have  a  (1  —  e)  optimal  solution  (if 
LB  >  (1  —  e)UB)  or  continue  to  choose  a  problem  with 
the  maximum  upper  bound  (Problem  3  in  Fig.  2(b))  and 
perform  partitioning  for  this  problem.  By  choosing  a 
problem  with  the  maximum  upper  bound  for  partition, 
we  can  ensure  that  UB  is  decreased  after  each  partition. 

An  important  technique  in  branch-and-bound  is  that 
we  can  remove  some  problems  from  further  consider¬ 
ation  and  thus  reduce  complexity.  In  particular,  if  we 
find  a  problem  z  with  LB  >  {I  —  e)UBz  (see  problem  4 
in  Fig.  2(c)),  we  conclude  that  this  problem  can  be  re¬ 
moved  from  further  consideration  without  loss  of  (1  — e) 
optimality. 

Figure  3  shows  the  main  algorithm.  Since  our  core 
optimization  space  is  finite  (with  finite  number  of  core 
variables  x")  and  gj",  and  each  core  variable  has  a  finite 
integer  value  set),  the  branch-and-bound  algorithm  is 
guaranteed  to  converge  (even  for  e  =  0)  [25].  Here, 
branch-and-bound  is  much  faster  than  brute-force  ex¬ 
haustive  search  because  non-improving  problems  are 
being  removed  during  the  process  to  avoid  wasting 
precious  cycle  time  in  future  computation.  As  a  result, 
for  all  the  network  instances  studied  in  Section  8,  brute- 
force  exhaustive  search  carmot  find  an  optimal  solution, 
while  our  algorithm  can  find  near-optimal  solutions  for 
all  network  instances  (with  various  network  sizes  and 
user  sessions). 

6  Determining  Bounds  and  Partitioning 
Problems 

The  main  algorithm  was  presented  in  Fig.  3.  Several 
components  (i.e.,  determining  upper  and  lower  bounds, 
partitioning  problems)  in  this  main  algorithm  are  yet  to 
be  developed.  These  algorithms  should  exploit  problem 
specific  properties  to  optimize  performance.  In  this  sec¬ 
tion,  we  design  algorithms  for  these  components. 

6.1  Determining  Upper  Bound 

To  find  an  upper  bound  for  a  problem  in  branch-and- 
bound  (see  lines  4  and  14  in  Fig.  3),  we  propose  to 
construct  a  linear  relaxation.  That  is  to  say,  we  linearize 


7 


Max  K 

l^keJ\r,k^i  -^ki  “T  J- 

(i  G  Mi) 

-  Qx^  <  0 

{i,j  eN,i^  j,me  Mij) 

1 

(i  e  Af,  m  e  Mi) 

m  1  oi.  ■n’"  0-  o’"  —  0 

{i,j  e  M,i  ^  j,m  Cl  Mij) 

Linear  constraints  for  ufjf. 

(i,j,  k  C  Af,i  j,m  G  Mij,m  S  Mk) 

<  0 

ij  ij  — 

{i,j  cN,ij^j,mC  Mij) 

(lCC,i  =  s{l)) 

“  J2ke)i\k^i,d(i)  /fehO  =  0 

{lCC,iCjV,ij^  s{l),d{l)) 

S^s(l)j^j,d(l)^i  ,  ...  vv  ^  r, 

l^iec  •'  In  2  Si  ^  ^ 

{i,j  cN,i^  3,  Mij  0) 

Linear  constraints  for  cYi 

(iJ  C  M,i  ijt  j,m  C  Mij) 

1  ^  ij  '  ij  '  ijk  — 

(i,j,k  cN,i  i=-  j,m  e  Mij,m  S  Mk) 

KJifl)  >  0 

(/  e  C,i,j  d{l),j  jh  s{l),Mij  0) 

(x,  q)  e  Q.Z 

Fig.  5.  Linear  reiaxation. 


Fig.  4.  A  convex  huii  for  =  ln(l  +  s™). 

all  constraints  in  Fig.  1  so  that  the  relaxed  problem  can  be 
solved  by  an  LP.  This  solution  provides  an  upper  bound. 

Note  that  in  Fig.  1,  and  log2(l+s^)  are  nonlinear 
terms.  For  product  we  apply  a  novel  method 

based  on  Reformulation-Linearization  Technique  (RLT)  [25]. 
This  is  done  by  introducing  a  new  variable  =  t^sfj 
and  add  linear  constraints  for  the  new  variable.  Suppose 
and  s™  are  bounded  by  (t™)L  <  and 

{sfj)L  <  sfij  <  {sfj)u,  respectively.  Thus,  we  have 

W  -  {tk)L\  ■  K  -  ^  0  ’ 

[{tk)u  -  C]  ■  Wj  -  (5™)l]  >  0  , 

-  C]  ■  [(s™)r/  -  s™]  >  0  . 

Substituting  =  t'f's"-,  we  have  the  following  linear 
consfrainfs  for 

{tT)L  ■  +  (s™)z,  •  C  -  <  (C)i  •  . 


y^k  )u  *  -^ij  +  -  U--J.  >  yt^  )u  •  ys^^)L  , 

y^k  )l  *  +  ysijju  ■  Cfc  -  U--J.  >  yt^  }l  •  ys^jju , 

fk  )u  ■  Sij  +  [Sij)u  ■  Cfe  -  Wijfc  <  (Cfe  )[/  •  [Sij)u  ■ 

For  fhe  nonlinear  ferm  log2(l+s[")  =  ^  ln(l+s[^),  we 

propose  fo  employ  fhree  fangenfial  supporfs  for  ln(l  + 
sfj),  which  is  a  convex  hull  linear  relaxafion  (see  Fig.  4). 
Suppose  fhaf  we  have  fhe  bounds  for  5™,  i.e.,  < 

sfj  <  We  introduce  a  variable  =  ln(l  +  s[^)  and 

consider  how  to  get  a  linear  relaxation  for  c™.  The  curve 
of  cfj  =  ln(l+s[j)  can  be  bounded  by  four  segmenfs  (or  a 
convex  hull),  where  segmenfs  I,  II,  and  III  are  fangenfial 
supporfs  and  segmenf  IV  is  fhe  chord  (see  Fig.  4).  In 
parficular,  fhe  fhree  fangent  segmenfs  are  fangenfial  at 
points  (1  +  (s[^)L,ln(l  +  (s[")l)),  (1  + /3, ln(l  + /?)),  and 
(1  +  (s™)c/,  ln(l  +  (s^)i/)),  where 

[l  +  (5-)^].[l  +  (4)^].[ln(l  +  (3[-)c,)-ln(l  +  (3-)^)]  ^ 

(s7])u  -  (s^)l 

is  the  horizontal  location  for  fhe  poinf  fhaf  is  infersecfed 
by  exfending  segmenfs  I  and  III;  segmenf  IV  is  fhe 
segmenf  fhaf  joins  poinfs  (1  +  (s™)L,ln(l  +  {s^)l))  and 
(l  +  (s'”)[/,ln(l  +  (s”- )[/)).  The  convex  region  defined  by 
fhe  four  segmenfs  can  be  described  by  fhe  following  four 
linear  consfrainfs. 

ll  +  is7])L]  ■  cT,-sT^<[l  +  {sTi)L]  ■  [ln(l  +  (s-)i)-l]  +  l  , 

{1 -\-  P)  ■  df  —  sfj  <  {1  -\-  fi)  ■  [ln(l  +  /3)  —  1]  +  1  , 

[l  +  (s]i)c/]  •  c^-s7j<[l  +  {s^)u]  ■  [ln(l  +  (s™)a)-l]  +  l  , 
[{s7])u-is^)L]  ■  cTi  +  [ln(l  +  (s™)L)-ln(l  +  (s™)c;)]  •  sTj 
>  {sTi)u  ■  ln(l  +  {sTj)L)  -  {sTj)L  ■  ln(l  +  {sTj)u)  ■ 

As  a  resulf,  fhe  nonlinear  In  (or  log)  ferm  is  relaxed  info 
linear  consfrainfs. 


After  relaxing  all  nonlinear  terms  for  a  problem,  say 
problem  z,  we  have  a  relaxed  problem  z  in  Fig.  5,  which 
is  an  LP.  In  problem  z,  x  and  q  are  the  vectors  for  all 
and  9™  variables,  respectively,  {q^)L, 

and  {q^)u  are  constant  bounds,  and  =  {(x,q)  : 
{x^)l  <  X™  <  {x^)u,{q?})L  <  C  <  (q^)u}  is  the  core 
optimization  space  of  (x,  q). 

A  relaxed  problem  z  can  be  solved  by  an  LP  in 
pol5momial  time.  Denote  its  solution  as  LP{z).  This  gives 
us  an  upper  bound  to  problem  z. 

6.2  Determining  Lower  Bound 

To  find  a  lower  bound  for  problem  z,  it  is  sufficient  to 
find  a  feasible  solution  to  this  problem.  By  feasible  solu¬ 
tion,  we  mean  that  it  satisfies  all  constraints  for  problem 
z  in  Fig.  1,  despite  that  the  objective  value  corresponding 
to  the  feasible  solution  may  not  be  optimal  (maximum). 
Although  any  feasible  solution  to  problem  z  can  serve 
as  a  lower  bound,  we  strive  to  find  one  feasible  solution 
that  can  offer  a  tight  lower  bound.  We  can  find  such 
a  feasible  solution  (denoted  as  ijjz)  based  on  LP{z)  by 
searching  its  neighborhood,  which  we  call  local  search. 

The  local  search  algorithm  begins  with  an  initial  fea¬ 
sible  solution.  Such  a  solution  may  be  far  away  from 
the  optimum  and  may  not  provide  a  tight  lower  bound. 
Thus,  we  will  iteratively  improve  the  current  solution  to 
achieve  a  better  lower  bound.  Until  we  can  no  longer 
improve  (increase)  the  lower  bound,  we  are  done. 

To  obtain  an  initial  feasible  solution,  we  set  x™  = 
(x™) L  for  scheduling  and  g™  =  (g™) l  for  power  control. 
Then  we  can  compute  SINK  value  by  (3).  When  an 
SINK  value  is  larger  than  or  equal  to  a,  the  achieved 
capacity  is  lUlog2(l  -I-  s™).  Otherwise  (i.e.,  SINK  <  a), 
the  transmission  is  considered  unsuccessful.  Note  that 
although  the  flow  rates  fij{l)  in  the  relaxed  solution 
LP{z)  guarantee  flow  balance  at  each  node,  such  flow 
rates  may  exceed  the  capacities  on  some  links  under 
the  initial  x("  and  oj"  values.  To  find  feasible  flow  rates 
under  current  x™  and  q^,  we  compare  the  achievable 
link  capacity  (under  current  xYj  and  g™  values)  to  the 
aggregated  flow  rates  fij{l)  on  each  link  z  — >  j  by 
computing  the  ratio  between  the  two  (denoted  as  A^) 
as  follows. 

_E^C>t.,W^l0g2(l  +  ^g) 

If  Xij  <  1  for  some  link  i  j,  then  the  aggregated 
flow  rates  exceed  the  link  capacity  and  the  link  capacity 
constraint  on  z  — >  j  is  violated.  In  this  case,  we  need  to 
scale  down  the  flow  rates  on  link  z  — >  j  (to  satisfy  link 
capacity  constraint)  and  the  flow  rates  on  all  other  links 
(to  maintain  flow  balance  in  the  network)  by  a  value 
A  <  Xij.  On  the  other  hand,  we  want  to  have  a  A  as 
large  as  possible  so  as  to  maximize  the  scaling  factor 
(our  objective).  Such  a  value  is  the  bottleneck  value  Xij 
among  all  links  (denoted  as  Amin  =  min{Ay  :  i,j  G 
Af,i  ^  j,Mij  ^  0}).  We  now  have  a  complete  solution 
'^min  '  /ij  (l),  (x”j')L,  (<?,7)i  foi'  routing,  scheduling,  and 
power  control,  respectively.  The  achieved  objective  is 


Initialization: 

1.  Set  xY]  =  (x"!j)L  and  =  {qY])L- 

2.  Compute  the  ratio  Xij  by  (17)  for  each  link  i  ^  j  and 
denote  Amin  =  min{Aij  :  i,j  &  J\f,i  j,  Mij  yf  0}. 

Iteration: 

3.  Suppose  Xij  =  Amin- 

4.  If  we  can  increase  q"j  on  a  used  band  { 

5.  Suppose  band  m  has  the  largest  q™  values  in  solution 
LP{z)  among  these  bands. 

6.  Increase  q^^  such  that  q™  <  (qY^)u  and  for  any  other 

7.  link  k  ^  h,  their  newly  updated  Xkh  >  Amin-  } 

8.  else,  if  we  can  increase  q^j  on  an  available  and  unused 
band  { 

9.  Suppose  band  m  has  the  largest  q™  values  in  solution 
LP{z)  among  these  bands. 

10.  Increase  qV}  such  that  qV}  <  {qYj)u  and  for  any  other 

11.  link  k  ^  h,  their  newly  updated  X^h  >  Amin- 

12.  If  qTt  increases,  then  xV^  =  1-1 

13.  else  the  iteration  terminates. 

Fig.  6.  Pseudocode  of  proposed  local  search  algorithm. 

Amin  •  at,  where  K  is  the  objective  value  in  the  relaxed 
solution  LP{z). 

In  the  next  iteration,  we  aim  to  improve  the  current 
solution.  Note  that  if  we  can  increase  Amin/  then  the 
current  solution  is  improved.  Suppose  link  z  — >  j  is  the 
link  with  A^  =  Amin-  To  increase  A^,  we  try  to  increase 
transmission  power  q™  on  some  band  m  under  the 
constraint  q”j  <  {q™)u-  Based  on  the  constraints  in  Fig.  1, 
we  may  update  the  values  of  other  variables  to  maintain 
feasibility.  For  example,  by  the  first  constraint  in  Fig.  1, 
we  need  to  increase  xj"  from  0  to  1  if  o’”  is  increased 
form  0  to  a  positive  value.  Moreover,  as  a  consequence 
of  increased  q”-,  the  interference  to  other  transmissions 
on  band  m  is  increased  and  thus  the  achieved  capacities 
for  other  links  are  decreased.  Thus,  gj”  can  be  success¬ 
fully  increased  only  if  for  any  other  link  fc  — >  h,  its 
updated  Xkh  will  not  fall  below  the  current  Amin-  h  the 
current  solution  can  be  improved  (with  a  larger  Amin)/ 
then  we  continue  to  the  next  iteration  of  improvement. 
Otherwise,  the  local  search  algorithm  terminates.  The 
pseudocode  of  our  local  search  algorithm  is  given  in 
Fig.  6. 

6.3  Partitioning  Probiem 

Our  proposed  partitioning  approach  differs  from  that 
in  standard  branch-and-bound  procedure.  In  standard 
branch-and-bound  procedure,  partitioning  (see  line  12 
in  Fig.  3)  is  done  by  choosing  a  variable  with  the  largest 
relaxation  error  and  uses  its  value  in  the  relaxed  solution 
LP{z)  to  divide  its  value  set  into  two  smaller  sets.  The 
reason  of  this  approach  (with  the  largest  relaxation  error) 
is  that  such  a  variable  is  likely  to  lead  to  a  larger  gap 
between  upper  and  lower  bounds.  Thus,  we  should 
partition  its  value  set  such  that  the  relaxation  error  will 
became  smaller.  This  division  (on  value  set)  also  divides 
the  optimization  space  for  problem  z  into  two  smaller 
spaces,  which  result  in  two  new  problems  zi  and  Z2, 
respectively. 

Such  standard  partitioning  technique,  however,  does 
not  explore  any  problem  specific  property  on  choosing 
partition  variables.  We  find  that  if  we  weigh  the  sig¬ 
nificance  of  each  variable  when  choosing  a  partitioning 
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variable,  the  complexity  of  the  overall  algorithm  can  be 
decreased  significantly.  For  our  problem,  x  variables  are 
more  important  than  q  variables.  That  is,  we  should  first 
decide  whether  or  not  a  band  is  used  for  fransmission. 
Only  if  a  band  is  used,  fhen  we  furfher  decide  fhe 
fransmission  power  on  fhis  band.  Thus,  we  firsf  choose 
one  X  variable  for  partitioning  and  affer  x  variables  are 
all  defermined  (eifher  0  or  1),  we  consider  q  variables 
for  parfifioning. 

When  choosing  a  specific  xfj  variable  for  parfifioning, 
we  pick  fhe  one  wifh  fhe  maximum  relaxation  error, 
which  is  defined  as  min{a;^,  1  —  x™},  where  is  fhe 
relaxed  solution  in  LP{z)  and  is  nof  necessarily  an 
infeger.  Once  partitioned,  fhe  value  sef  for  fhis  x  can 
have  elemenf  0  or  1.  This  variable  can  be  replaced  by  a 
consfanf  0  or  1  in  fhe  two  new  subproblems. 

We  furfher  observe  fhaf  by  fixing  fhe  value  of  x™, 
some  ofher  variables  may  also  be  fixed.  We  nofe  fhaf 
by  (1),  a  node  carmof  fransmif  and  receive  on  fhe  same 
band.  Based  on  fhis  observafion,  if  fhe  value  of  x™  is  sef 
fo  1,  fhen  we  have  x^  =  0  for  k  G  Af,k  i,m  €  Mk  and 
x^  =  0  for  h  G  Af,h  ^  j,m  G  Mh-  On  fhe  ofher  hand,  if 
fhe  value  of  x™  is  sef  fo  0,  fhen  we  have  g*”  =  0  based 
on  (2). 

Affer  all  x  variables  are  fixed,  we  can  sfarf  parfifion 
on  q  variables.  Again,  we  choose  one  q  variable  wifh 
fhe  maximum  relaxafion  error  for  parfifioning.  The  re¬ 
laxation  error  of  a  variable  g™  is  defined  as  minlg™  — 
[g™J,  [g™J  -I-  1  —  g™},  where  g^  is  fhe  relaxed  solution 
in  LP{z)  and  is  nof  necessarily  infeger.  The  value  sef 
of  q^  in  problem  0  is  {(g’”)^,  (g’”)L  + 

Then  ifs  new  value  sef  in  fhe  fwo  subproblems  will 
be  {(g-)i,  (g™)L  +  1,  •  •  • ,  [g™]}  and  {[g^J  +  1,  Lg“J  + 
2,  •  •  • ,  (g,“)c/}/  respectively 


7  Discussion 

We  now  discuss  how  fo  inferpref  and  apply  our  fheoref- 
ical  resulfs.  Nofe  fhaf  in  previous  secfions,  fhe  problem 
setting  under  sfudy  and  fhe  proposed  solution  procedure 
are  concerned  wifh  an  insfance  of  CRN  where  fhe  avail¬ 
able  bands  af  each  node  is  given  (sfafic).  Our  fheorefical 
resulf  offers  a  performance  benchmark  for  a  mulfi-hop 
CRN.  Such  a  resulf  is  nof  yef  available  in  fhe  liferafure. 

We  make  no  claim  fhaf  our  solution  will  be  deployed 
in  ifs  currenf  form  in  fhe  real  world.  Insfead,  we  po¬ 
sition  our  resulf  fo  be  used  as  a  reference  opfimal  re¬ 
sulf /performance  benchmark  for  any  ofher  algorifhms 
fhaf  will  be  developed  for  deploymenf.  To  sfudy  fhe 
performance  of  ofher  proposed  algorifhm,  we  can  record 
ifs  performance  over  a  time  period  T.  For  performance 
benchmarking/comparison,  we  can  run  our  solufion 
offline  by  breaking  up  fhe  time  period  T  info  smaller 
time  infervals  where  fhe  available  bands  can  be  assumed 
sfafic  wifhin  each  small  time  inferval  '  •  For  each 

time  inferval  ti,  our  algorifhm  can  provide  a  solufion  tt^. 
Then  a  complefe  solufion  If  for  fhe  entire  time  period 
T  is  obfained  by  piecing  up  fogefher  all  solutions  tt^. 
Solufion  If  is  eifher  opfimal  (by  setting  £  =  0)  or  near- 


opfimal  (if  e  >  0).^  Then  we  can  compare  how  far  fhe 
performance  of  fhe  ofher  proposed  algorifhm  is  from  fhe 
opfimum. 

8  Numerical  Results 

In  fhis  secfion,  we  presenf  numerical  resulfs  for  fhe 
proposed  solufion.  Our  goals  are  fo  demonsfrafe  fhe 
efficacy  of  fhe  solufion  procedure  and  fo  offer  quanfifa- 
five  undersfanding  on  fhe  joinf  opfimizafion  af  differenf 
layers  in  fhe  SINR  model. 

8.1  Simulation  Setting 

We  consider  20-,  30-,  and  50-node  CRNs  with  each  node 
randomly  located  in  a  50x50  area  (see  Tables  1,  3,  and 
8).  For  the  ease  of  exposition,  we  normalize  all  units  for 
distance,  bandwidth,  rate,  and  power  based  on  (8)  with 
appropriate  dimensions.  We  assume  there  are  \M  \  =  10, 
20,  and  30  frequency  bands  in  the  network  for  the  20-, 
30-,  and  50-node  CRNs  and  each  band  has  a  bandwidth 
of  W  =  50.  At  each  node,  only  a  subset  of  these  bands 
is  available.  For  the  20-  and  30-node  CRNs,  we  assume 
there  are  5  user  communication  sessions  (see  Tables  2 
and  4)  and  for  the  50-node  CRN,  the  number  of  user 
communication  sessions  is  10  (see  Table  9).  The  source 
node  and  destination  node  for  each  session  are  randomly 
selected  and  the  minimum  rate  requirement  of  each 
session  is  randomly  generated  within  [1, 10]- 

We  assume  =  d”"*  and  the  SINR  threshold  a  =  3 
[12].  We  assume  that  under  the  maximum  transmission 
power,  a  node  at  distance  20  can  receive  data  when  there 

is  no  interference.  Thus,  we  have  =  a,  i.e.,  the 

maximum  transmission  power  is  Pmax  =  o  •  (20)^  •  rjW  = 
4.8  •  lO^rjW.  We  assume  that  power  control  can  be  done 
in  Q  =  10  levels.  For  our  proposed  solution,  we  set  e  to 
0.1,  which  guarantees  that  the  solution  is  at  least  90% 
optimal. 

We  note  that  the  brute  force  approach  carmof  solve  the 
problem  even  for  the  20-node  CRN.  The  solution  space  of 
the  capacity  problem  in  Section  3.3  includes  all  possible 
sets  of  values  for  (x™,  g™,  K,  fij{l)).  Thus,  the  number  of 
solutions  examined  in  the  brute  force  approach  is  clearly 
more  than  the  number  of  all  possible  sets  of  values  for 
g™  variables.  We  now  analyze  the  number  of  all  possible 
sets  of  values  for  g™  variables  for  the  20-node  CRN.  The 
number  of  g™  variables  is  about  20  •  (20  —  1)  •  5  =  1900, 
where  5  is  an  approximation  of  the  average  number  of 
available  bands  on  a  link.  Each  g™  variable  has  (Q  -I- 1) 
or  11  possible  values.  Thus,  the  number  of  all  possible 
sets  of  values  for  g™  variables  is  about  11^®°°.  Therefore, 
for  the  20-node  CRN,  the  number  of  solutions  examined 
in  the  brute  force  approach  is  at  least  11^®°°.  Suppose 

1.  We  discuss  the  case  of  e  =  0.  The  discussion  for  e  >  0  is  similar 
and  thus  is  omitted.  First  of  all,  since  available  bands  are  static  in  any 
time  interval  ti,  solution  ni  is  optimal  as  we  discussed  in  Section  5.2. 
That  is,  there  is  no  solution  if  ^  that  can  achieve  better  performance  than 
TTi  in  time  interval  ti.  Since  this  claim  holds  for  each  time  interval  ti, 
there  is  no  solution  n  that  can  achieve  better  performance  than  IT  for 
any  time  interval.  Therefore,  solution  IT  is  optimal. 


TABLE  1 

Location  and  available  frequency  bands  at  each  node  for  a  20-node  network. 
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Node 

Location 

Available  Bands 

Node 

Location 

Available  Bands 

1 

(0.1,  9.9) 

1,  2,  3,  4,  7,  8,  9,  10 

11 

(28.1,25.6) 

1,  2,  3,  4,  5,  6,  7,  8,  9,  10 

2 

(29.2,31.7) 

1,  2,  3,  4,  5,  7,  8,  10 

12 

(32.3,38) 

1,  8,  9,  10 

3 

(3,31.1) 

1,  4,  5,  6 

13 

(47.2,2.6) 

3,  5,  10 

4 

(11.8,40.1) 

1,  2,  3,  4,  6,  9,  10 

14 

(44.7,15) 

2,  3,  6,  7,  8 

5 

(15.8,9.7) 

1,  2,  3,  5,  6,  8,  9 

15 

(44.7,24) 

1,  2,  3,  4,  5,  6,  7,  8,  9,  10 

6 

(16.3N9.5) 

3,  5,  6,  8,  9 

16 

(47.9N3.8) 

1,  3 

7 

(0.6,27.4) 

1,  4,  8,  9,  10 

17 

(46.4N6.8) 

1,  7,  9 

8 

(22.6,40.9) 

1,  2,  3,  5,  7,  9,  10 

18 

(11.5,12.2) 

2,  5,  6,  10 

9 

(35.3,10.3) 

2,  9 

19 

(28.2,14.8) 

4,  5,  6,  7,  8,  9,  10 

10 

(31.9,19.6) 

1,  2,  3,  4,  5,  6,  7,  8,  9,  10 

20 

(2.5,14.5) 

1,  7,  10 

TABLE  2 

Source  node,  destination  node,  and  minimum  rate 
requirement  of  each  session  in  the  20-node  network. 


Session 

1 

Source  Node 
s{l) 

Dest.  Node 
d{l) 

Min.  Rate  Req. 
r{l) 

1 

16 

10 

9 

2 

18 

3 

1 

3 

12 

11 

4 

4 

13 

17 

3 

5 

15 

6 

2 

50 
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Fig.  7.  The  routing  topology  for  the  20-node  5-session 
network. 

one  solution  can  be  examined  in  10“®  second.  Then  the 
running  time  is  ^  2^2^1894  seconds,  which 

is  111®94/(305  .  24  •  60  •  60)  >  llis94/208  >  yg^rs. 

Therefore,  the  brute  force  approach  carmof  be  used  fo 
solve  our  problem  even  for  the  20-node  CRN. 

8.2  Results 

For  the  20-node  network  with  5  sessions,  our  so¬ 
lution  achieves  a  scaling  factor  K  =  13.24.  Then, 
based  on  the  minimum  rate  requirement  r{l)  in  Ta¬ 
ble  2,  the  flow  rates  K  ■  r{l)  for  the  five  sessions 
are  119.16, 13.24, 52.96, 39.72, 26.48,  respectively.  Figure  7 
shows  the  routing  topology  in  the  final  solution.  The 
flow  rates  for  each  session  on  the  links  along  its  path 


are  as  follows. 

Session  I  =  1:  /le  12(1)  =  119.16,  /i2  8(l)  =  103.30, 
/i2.ii(l)  =  15.86,  /8,2(1)  =  103.30,  /ii,io(l)  =  15.86, 
/2.io(l)  =  103.30; 

Session  I  =  2:  /i8.2o(2)  =  13.24,  /2o.i(2)  =  13.24, 
/i,7(2)  =  13.24,  /7,3(2)  =  13.24; 

Session  I  =  3:  /i2,ii(3)  =  52.96; 

Session  I  =  4:  /i3’i4(4)  =  39.72,  fu  17(4)  =  39.72; 

Session  I  =  5:  /i5U9(5)  =  26.48,  fig^eib)  =  26.48. 

It  is  easy  to  verify  that  flow  balance  holds  at  source 
and  destination  nodes  for  each  session,  as  well  as  at 
all  relay  nodes.  Note  that  flow  splitting  and  multi-path 
routing  are  used  for  session  1,  which  has  the  largest  rate 
requirement. 

Our  solution  also  solves  the  scheduling  variables  x™ 
as  follows.  Note  that  it  is  sufficient  to  list  only  non-zero 
x")  variables. 

‘‘J 

Band  to  =  1:  3  =  1,  x\q  32  =  li 

Band  to  =  2:  a;|  2  =  li 
Band  to  =  3:  a;f3  14  =  1; 

Band  to  =  4:  7  =  1,  X2  ig  =  1; 

Band  to  =  5:  40  =  1' 

Band  to  =  6:  19  =  1; 

Band  to  =  7:  3^14^17  =  1,  2:20,1  =  If 
Band  to  =  8:  a;i2,ii  =  1; 

Band  to  =  9:  a;i2,8  =  1,  a;fg  g  =  1; 

Band  to  =  10;  xjg  20  =  1- 

The  transmission  power  levels  on  their  respective  fre¬ 
quency  bands  are  as  follows. 

Band  to  =  1:  3  =  1,  gjg_i2  =  7; 

Band  to  =  2:  2  =  2; 

Band  to  =  3:  44  =  2; 

Band  to  =  4:  y  =  7,  10  =  2; 

Band  to  =  5:  49  =  1; 

Band  to  =  6:  (745,49  =  9; 

Band  to  =  7:  9[4_47  =  1,  g|o,i  =  If 
Band  to  =  8:  (742,41  =  3; 

Band  to  =  9:  g?2,8  =  1/  9i9,6  =  3f 
Band  to  =  10:  (713,20  =  1- 

Note  that  the  same  frequency  band  may  be  used  by 
concurrent  transmissions.  For  example,  since  Xy  3  =  1 
and  Xyg  42  =  1,  we  have  that  both  nodes  7  and  16  are 
transmitting  on  band  1.  Such  concurrent  transmissions 


TABLE  3 

Location  and  available  frequency  bands  at  each  node  for  a  30-node  network. 
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Node 

Location 

Available  Bands 

Node 

Location 

Available  Bands 

1 

(7,0.7) 

1,  2,  6,  7,  16,  17,  19,  20 

16 

(30.3,28.1) 

7,  8,  11,  16,  17,  19,  20 

2 

(5,4) 

3,  5,  9,  12,  14,  15 

17 

(32,41.1) 

7,  11,  16,  17,  19,  20 

3 

(6.8,14) 

1,  2,  6,  7,  8,  11,  16,  17,  19,  20 

18 

(14.1,33.7) 

3,  4,  5 

4 

(15.7,3.3) 

1,  2,  7,  16,  20 

19 

(23,46.4) 

3,  12,  15 

5 

(9.5,17) 

3,  4,  5,  9,  12 

20 

(30.3,9.3) 

5,  9 

6 

(19.4,17.1) 

1,  2,  6,  7,  8,  16,  19,  20 

21 

(17.6,29.2) 

1,  2,  6,  7,  8,  11,  16,  17,  19,  20 

7 

(34.7, 14.6) 

3,  4,  5,  9,  12,  14 

22 

(27.1,27.8) 

9,  12,  14,  15 

8 

(4.9,25.9) 

3,  4,  12 

23 

(26.9,45.9) 

3,  4,  5,  9,  10,  12,  13,  14,  15,  17 

9 

(46.6,42.1) 

10,  18 

24 

(43.3,32.4) 

1,  2,  11,  16,  17,  20 

10 

(8.3,38.3) 

3,  4,  5,  9,  14 

25 

(45.4,8.2) 

3,  4,  5,  9,  12,  14 

11 

(26.7,11.1) 

1,  6,  7,  8,  11,  16,  17,  19,  20 

26 

(43.4,35) 

3,  5,  9,  15 

12 

(36.4N7.3) 

10,  13,  18 

27 

(41.3N5.1) 

1,  16,  20 

13 

(24.3,21.2) 

1,  2,  6,  8,  11,  19 

28 

(14.4,30.3) 

1,  2,  6,  7,  8,  11,  16,  17,  20 

14 

(23.1,0.8) 

3,  5,  9,  14 

29 

(41.6,41.7) 

3,  4,  5,  9,  10,  12,  14,  15,  18 

15 

(21.4, 19.2) 

4,  9,  12,  14 

30 

(25.9,12) 

1,  2,  6,  7,  8,  11,  16,  17,  19,  20 

TABLE  4 

Source  node,  destination  node,  and  minimum  rate 
requirement  of  each  session  in  the  30-node  network. 


Session 

1 

Source  Node 
s{l) 

Dest.  Node 
d{l) 

Min.  Rate  Req. 
r{l) 

1 

16 

28 

4 

2 

24 

11 

7 

3 

13 

1 

1 

4 

19 

29 

8 

5 

26 

15 

1 

Fig.  8.  The  routing  topology  for  the  30-node  5-session 
network. 

are  allowed  as  long  as  the  SINR  at  each  receiving 
node  is  no  less  than  a.  For  example,  SINR  at  receiv¬ 
ing  node  12  on  band  1  is  sjg  12  =  = 

f’max+*7,i2''?7.3 

2.083-id-°+(8°0ii-i0-^)-i  =  4-22,  which  is  larger  than  a  =  3. 
Thus,  transmission  16  ^  12  on  band  1  is  successful. 
Following  the  same  token,  we  can  compute  all  SINR 
values  as  follows. 

Band  to  =  1:  ^  =  118.47,  =  4.22; 

Band  to  =  2:  s|  2  =  5.84; 


Band  to  =  3:  sfg  =  3.75; 

Band  to  =  4:  sf  7  =  3.33,  82,10  =  3.14; 

Band  to  =  5:  sfuo  =  18.87; 

Band  to  =  6:  sfg  79  =  3.39; 

Band  to  =  7:  S74  77  =  1261.14,  820,1  =  155.46; 

Band  to  =  8:  842,11  =  4.90; 

Band  to  =  9:  842,8  =  3.56,  qfg  Q  =  4.74; 

Band  to  =  10:  84°  20  =  6.45. 

As  expecfed,  we  see  fhaf  fhe  achieved  SINR  af  each 
receiving  node  on  each  band  is  larger  fhan  a  =  3  in 
our  solufion. 

For  each  link,  we  can  furfher  verify  fhaf  flow  rafes 
on  fhis  link  do  nof  exceed  ifs  capacify.  For  example,  for 
link  16  ^  12,  fhere  is  a  flow  rafe  /i6,i2(l)  =  119.16  on 
fhis  link.  The  achieved  capacify  is  FFlog2(l  -I-  849  42)  = 
FFlog2(l -L  4.22)  =  119.21. 

The  above  resulfs  are  for  20-node  nefwork.  The  resulfs 
for  30-  and  50-node  nefworks  are  similar  and  we  abbre- 
viafe  our  discussion.  For  fhe  30-node  5-session  nefwork, 
fhe  available  bands  af  each  node  and  fhe  locafion  of 
each  node  are  shown  in  Table  3.  The  source,  desfina- 
fion,  and  minimum  rafe  requiremenf  of  each  session 
are  shown  in  Table  4.  Our  solufion  gives  a  routing 
topology  in  Fig.  8.  The  achieved  scaling  factor  is  31.18. 
Then,  based  on  fhe  minimum  rafe  requiremenf  r(l)  in 
Table  4,  fhe  flow  rafes  K  ■  r{l)  for  fhe  five  sessions 
are  124.72,218.26,31.18,249.44,31.18,  respectively  The 
flow  rafes  for  each  session  on  links  along  ifs  paths  are 
shown  in  Table  5.  Our  solution  solves  the  scheduling 
variables  a;™  and  power  control  variables  g™.  Again,  we 
list  only  non-zero  x™  variables  in  Table  6  and  non-zero 
variables  in  Table  7. 

For  the  50-node  10-session  network,  the  available 
bands  at  each  node  and  the  location  of  each  node  are 
shown  in  Table  8.  The  source,  destination,  and  minimum 
rate  requirement  of  each  session  are  shown  in  Table  9. 
Our  solution  gives  a  routing  topology  in  Fig.  9.  The 
achieved  scaling  factor  is  13.36.  The  detailed  solution  on 
power  control,  scheduling,  and  routing  are  similar  to  the 
20-  and  30-node  networks  and  are  omitted  to  conserve 
space. 
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TABLE  5 

Flow  rates  for  each  session  in  the  30-node  5-session  network. 


Session  1 

Session  rate  K  ■  r{l) 

Flow  rates  of  session  1  on  each  link 

1 

124.72 

/i6,2i(l)  =  124.72,  /21.28(1)  =  124.72; 

2 

218.26 

/24.16(2)  =  113.91,  /24.i7(2)  =  104.36,  /i7,i6(2)  =  104.36,  /i6.2i(2)  =  24.22, 
/i6.13(2)  =  194.05,  /2i,28(2)  =  24.22,  /28,i3(2)  =  24.22,  /i3,ii(2)  =  160.59, 
/i3,3o(2)  =  57.68,  /3o.ii(2)  =  57.68; 

3 

31.18 

/i3,3o(3)  =  31.18,  /3o.4(3)  =  31.18,  /4,i(3)  =  31.18; 

4 

249.44 

/i9,23(4)  =  167.49,  /i9.29(4)  =  81.95,  /23.26(4)  =  81.26,  /23,29(4)  =  86.23, 
/26,29(4)  =  81.26; 

5 

31.18 

/26.22(5)  =  31.18,  /22.15(5)  =  31.18. 

TABLE  6 

Scheduling  for  the  30-node  network. 


Band 

Scheduling 

Band 

Scheduling 

Band 

Scheduling 

Band 

Scheduling 

1 

a^4,l  =  1/  ®21,28  = 

6 

®13,11  =  1/ 

11 

®17,16  =  1/ 

16 

*24,17  =  1/  *30,11  =  f/ 

2 

®28,13  = 

7 

®16,21  =  1/ 

12 

*19, 23  =  1/ 

17 

*24,16  =  1/ 

3 

®19,29  =  1/ 

8 

®16,13  =  1/ 

13 

*i,'  =  0; 

18 

*i,'  =  0; 

4 

®23,29  = 

9 

®26,22  =  1/ 

14 

*2ll5  =  1/ 

19 

*13,30  =  1/ 

5 

®26,29  =  1/ 

10 

H 

II 

O 

15 

*23,26  =  1/ 

20 

*30,4  —  I- 

TABLE  7 

Transmission  power  levels  for  the  30-node  network. 


Band 

Transmission  power 

Band 

Transmission  power 

Band 

Transmission  power 

Band 

Transmission  power 

1 

94,1  =  1/  921,28  =  1; 

6 

9i3,ii  —  2; 

11 

9i7,16  —  2; 

16 

924,17  —  3,  §30, 11  =  1; 

2 

^  Q 

928,13  —  'J/ 

7 

916,21  =  4; 

12 

919,23  =  1; 

17 

924,16  =  7; 

3 

919,29  =  9/ 

8 

916,13  =  2; 

13 

II 

o 

18 

o 

II 

2:? 

4 

923,29  =  4; 

9 

926,22  =  7; 

14 

922,15  —  1/ 

19 

913,30  —  1; 

5 

926,29  =  1/ 

10 

^:c 

II 

o 

15 

923,26  =  10; 

20 

930,4  =  4. 

TABLE  8 

Location  and  available  frequency  bands  at  each  node  for  a  50-node  network. 


Node 

Location 

Available  Bands 

Node 

Location 

Available  Bands 

1 

(11.1,21.7) 

2,  3,  4,  8,  25 

26 

(25.2,27.2) 

10,  14,  20,  24,  26 

2 

(0.1,4) 

6,  7,  10,  13,  14,  20,  23,  24,  26,  28 

27 

(22.5,42.2) 

5,  9,  12,  16,  18,  27,  29,  30 

3 

(7.2,16.6) 

6,  10,  14,  20,  23,  24,  26 

28 

(30,31.5) 

6,  13,  24,  26,  28 

4 

(11,32.2) 

6,  7,  10,  13,  14,  20,  23,  24,  26,  28 

29 

(35,22.1) 

6,  10 

5 

(16.3,3.6) 

10,  13,  14,  20,  23 

30 

(25.7,6.2) 

5,  9,  12,  17,  18,  22,  27,  29,  30 

6 

(14.5,24.7) 

8,  11,  25 

31 

(34.1,12.4) 

9,  12,  16,  17,  30 

7 

(14.9,13.7) 

5,  9,  12,  16,  17,  18,  22,  27,  29,  30 

32 

(26.4,30) 

5,  9,  12,  16,  17,  18,  22,  27,  29,  30 

8 

(19.5,14.9) 

7,  24,  28 

33 

(14.1,40.7) 

1,  2,  25 

9 

(26.6, 13.4) 

1,  19,  21,  25 

34 

(34.4,46.5) 

9,  17,  18,  30 

10 

(22.5,29.3) 

1,  3,  4,  8,  11,  15,  19 

35 

(19,22.5) 

1,  6,  7,  10,  13,  14,  20,  23,  24,  28 

11 

(24.6,40.5) 

3,  8,  25 

36 

(39.9,25.1) 

6,  13,  14,  20,  23,  24,  26,  28 

12 

(38.4,13.1) 

2,  8,  11,  15 

37 

(20.3,18.2) 

1,  2,  3,  4,  8,  11,  15,  19,  21,  27 

13 

(4,3.9) 

9,  12,  16,  22,  27,  29,  30 

38 

(10,20.5) 

6,  7,  10,  13,  14,  20,  23,  24,  26,  28 

14 

(6.1,18.6) 

9,  12,  16,  17,  18,  22,  27,  30 

39 

(20.5,21.4) 

1,  2,  3,  4,  8,  11,  15,  19,  21,  25 

15 

(38.5,22.6) 

2,  4,  11,  15,  19,  21,  25 

40 

(37.1,28.6) 

7,  10,  13,  14,  20,  23,  24,  26 

16 

(1.2,24.3) 

5,  9,  12,  17,  22,  29,  30 

41 

(44.1,16.1) 

1,  15,  21 

17 

(4.9,42.3) 

5,  27 

42 

(41.1,6) 

9,  29 

18 

(18.5,1.4) 

5,  9,  12,  17,  18,  27,  30 

43 

(43, 18.8) 

5,  9,  12,  16,  18,  22 

19 

(16.9,29.1) 

3,  4,  10,  11,  12,  15 

44 

(45.4,24.2) 

9,  12,  16,  17,  18,  30 

20 

(33.5, 10.4) 

7,  13,  14,  20,  23,  24,  26,  28 

45 

(36.2,41.2) 

5,  9,  17,  27,  29,  30 

21 

(25.6,12.8) 

6,  7,  20,  23,  24,  28 

46 

(27.5,32.3) 

12,  16,  17,  18,  29,  30 

22 

(45.2,45.5) 

2,  8,  15,  19 

47 

(47.8,13.8) 

22,  27,  29,  30 

23 

(43.6,22.7) 

1,  2,  3,  4,  11,  15,  19,  21 

48 

(8.9,14.8) 

5,  30 

24 

(10.6,40.5) 

4,  15,  19,  21,  25 

49 

(6.8,  6.2) 

5,  9,  12,  16,  17,  27,  30 

25 

(18.2,32.7) 

9,  12,  18,  22,  27 

50 

(11.7,35.8) 

1,  2,  3,  4,  8,  11,  15,  19,  21,  25 

13 


TABLE  9 

Source  node,  destination  node,  and  minimum  rate 
requirement  of  each  session  in  the  50-node  network. 


Session 

1 

Source  Node 
s(Z) 

Dest.  Node 
d(l) 

Min.  Rate  Req. 

r{l) 

1 

21 

4 

4 

2 

5 

26 

7 

3 

19 

20 

6 

4 

33 

6 

10 

5 

37 

10 

9 

6 

23 

11 

2 

7 

25 

46 

3 

8 

42 

43 

9 

9 

44 

27 

8 

10 

47 

30 

1 

50 
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Fig.  9.  The  routing  topoiogy  for  the  50-node  10-session 
network. 

In  summary,  we  have  the  following  observations  for 
our  numerical  results.  To  maximize  capacity,  it  is  neces¬ 
sary  to  consider  multi-path  routing  (i.e.,  to  allow  flow 
splitting)  and  to  consider  spatial  reuse  (i.e.,  to  allow 
concurrent  transmissions  to  use  the  same  band).  Under 
multi-path  multi-hop  routing,  even  the  number  of  user 
sessions  are  small,  the  number  of  nodes  in  the  network 
that  will  be  affect  by  these  sessions  (i.e.,  multi-hop 
relays)  can  be  large. 

9  Conclusion 

In  this  paper,  we  investigated  a  capacity  problem  for 
a  multi-hop  cognitive  radio  network  under  the  SINR 
model.  We  showed  the  complex  inter-relationship  among 
power  control,  scheduling,  and  routing  in  this  problem 
and  developed  a  mathematical  formulation  for  joint 
optimization.  We  presented  a  solution  to  this  problem 
based  on  branch-and-bound  framework.  The  novelty  of 
our  contribution  is  the  determination  of  the  core  opti¬ 
mization  space  for  this  problem,  which  is  much  smaller 
than  the  whole  optimization  space,  and  the  design  of 
several  efficient  algorithms  within  branch-and-bound 
framework  that  exploit  the  physical  significance  of  core 
variables  to  obtain  tighter  upper  and  lower  bounds.  The 


proposed  solution  is  able  to  guarantee  (1  —  e)  optimal 
solution  and  can  be  used  as  a  performance  benchmark 
for  other  proposed  algorithms  in  the  future. 
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